Quantum Dot as a Spin— Current Diode: a spin-dependent transport study 
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We report a study of spin dependent transport in a system composed of a quantum dot coupled to 
a normal metal lead and a ferromagnetic lead (NM-QD-FM). We use the master equation approach 
to calculate the spin-resolved currents in the presence of an external bias and an intra-dot Coulomb 
interaction. We find that for a range of positive external biases (current flow from the normal metal 
to the ferromagnet) the current polarization p — (If— ij.)/ (It+I±) is suppressed to zero, while for the 
corresponding negative biases (current flow from the ferromagnet to the normal metal) p attains a 
relative maximum value. The system thus operates as a rectifier for spin-current polarization. This 
effect follows from an interplay between Coulomb interaction and nonequilibrium spin accumulation 
in the dot. In the parameter range considered, we also show that the above results can be obtained 
via nonequilibrium Green functions within a Hartree-Fock type approximation. 

PACS numbers: PACS number 



I. INTRODUCTION 

Polarized transport in spin-dependent nanostruc- 
tures is a subject of intense study in the emerg- 
ing field of spintronicsji due to its relevance to 
the development of novel spin-based devices £££ In 
addition, transport through QDs provides informa- 
tion about fundamental physical phenomena in spin- 
dependent and strongly correlated systems, such as the 
Kondo effect£&£&2 the Coulomb- and spin-blockade 
effects ) 10 ' 11 ' 12 ' 13 ' 14 spin valve effect and tunnelling mag- 
netoresistance (TMR)^^L 18 ' 19 ' 20 ' 21 ' 22 ' 23 ' 24 ' 25 ' 26 etc. 
Novel spin filters and pumps have also been proposed 
using QDs coupled to normal metal leads. 27 ! 28 ' 29 A sys- 
tem of particular interest in this context comprises a 
quantum dot or a metallic grain coupled to ferromag- 
netic leads. The ferromagnetism of the leads intro- 
duces spin-dependent tunneling rates between the leads 
and the central region. This results in a nonzero net 
spin in the central region for asymmetric magnetiza- 
tion geometries. This effect is called spin accumula- 
tion or spin imbalance. 30 ' 31 ' 32 It has been shown that 
spin accumulation affects several transport properties, 
such as magnetoresistance j 18 ' 35 (negative) differential 
resistanc o 10 ' 35 and the zero-bias anomaly^ ' 11 ' 18 In addi- 
tion it provides a way to generate and control the current 
spin polarization via gates or bias voltages ; 33 ' 34 which is 
one of the main goals within spintronics. 

Systems composed of a nonmagnetic lead and a ferro- 
magnetic lead with a quantum dot or a quantum wire 
as spacer have been analyzed recently. It was pointed 
out that if the spacer is a dot and the ferromagnetic lead 
is half-metallic, a novel mesoscopic current-diode effect 



arises 



4.35.36 



Spin-current rectification was also predicted 
in an asymmetric system composed of a ferromagnetic 



(Fc or Ni) and nonmagnetic (Au or Pd) contacts cou- 
pled to each other via a molecular wire^ 7 - Additionally, 
it was pointed out that a NM-QD-FM system can op- 
erate as a spin-filter and as a spin-diode^ 8 - In Ref. [38[ 
the authors use the bias voltage to change the resonance 
position of the dot level with respect to the spin-split 
density of states of the ferromagnetic lead. This gives 
rise to spin-dependent currents. 

Here we study spin-resolved currents in a single-level 
quantum dot attached to a nonmagnetic lead ( "left lead" ) 
and to a ferromagnetic lead ("right lead"), Fig. 1. As 
we shall show, the magnetic asymmetry between the left 
and right terminals results in a rectification of the current 
polarization for a particular bias range for which the sin- 
gle electron channel is on-resonance, and the double- 
occupancy channel + U is off-resonance. More pre- 
cisely, when the nonmagnetic lead operates as an emitter 
and the ferromagnetic lead as a collector, defined as the 
positive bias (eV > 0), the current is unpolarizcd. In 
contrast, when the ferromagnetic lead is the emitter and 
the nonmagnetic lead is the collector (negative bias) a 
spin-polarized current arises. Importantly, this rectifica- 
tion occurs only in this particular bias range, as we shall 
demonstrate both analytically and numerically. This is 
attributed to an interplay between nonequilibrium spin 
accumulation and Coulomb interaction within the dot. 
For high enough bias voltages, the current polarization 
is essentially symmetric with respect to the bias, and no 
rectification is found. 

In the main body of the text we employ the master- 
equation approach of Glazman and Matvee\*2 9 . to describe 
the spin-dependent transport through the NM-QD-FM 
junction in the sequential tunneling regime (To -C fc^T,— 
where Tq is a characteristic tunneling rate). An alterna- 
tive description in terms of nonequilibrium Green func- 
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FIG. 1: System studied: a nonmagnetic quantum dot coupled 
via tunneling barriers to a nonmagnetic left lead and a ferro- 
magnetic right lead. A bias voltage V is applied across the 
system so that the left (^ll) and right (hr) chemical potentials 
differ by \il — ^r = eV . 



tions is also presented in the appendix, that corroborates 
our results obtained via master equation. 



II. MODEL AND MASTER EQUATION 
APPROACH 



The NM-QD-FM system we study is schematically il- 
lustrated in Fig. [T] An external bias voltage V drives the 
system away from equilibrium thus imposing a chemical 
potential imbalance between the left (L) and the right 
(R) leads: fiL — Pr = eV, where fiuR) is the chemical 
potential of the lead L(R) and e is the absolute value of 
the electron charge (e > 0). The system Hamiltonian is 



(1) 



kcrr; 

+ ^e d dld a + Ud^d^ 

a 
kcr?7 

where eko-?; is the free-electron energy with wave vector k 
and spin a in lead r\ (77 = L, R), e d is the spin-degenerate 
dot level, U is the on-site Coulomb repulsion and the 
operators Ck ari (c^. ) and d a (eft) destroy (create) an 
electron with spin a in the lead 77 and in the dot, re- 
spectively. The matrix element t^^ gives the lead-dot 
coupling. We do not consider any spin-flip processes. 

To calculate the current we use rate equationS ) 24 i 39 
which yield 

where we have assumed H = 1. The parameter rJJ lcr 
corresponds to the rate of adding one electron to the 
dot coming from lead 77, and r^ 0(T is the rate of moving 
one electron from the dot to lead 77. In addition, Tq 1o . 
and r^ 0(T give the rates of moving one electron with spin 
a to and from the dot, respectively, when it is already 
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occupied by one electron with opposite spin. Following 
Ref. [24| we define n a = (h a ) and ri|j_ = (n-|-nj,) (n a = 
d^da) as the dot single and double average occupancies, 
respectively. The tunneling rates are 
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(4) 
(5) 
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where f v = l/{exp[(e d - ^/{ksT)] + 1} and f n = 
l/{exp[(e d + C/-/j., ) )/(fc B T)] + l}. The rates T% and rjj 
are related to the spin-resolved density of states of lead 
77 via T r < = 2n\t\ 2 Pail (e d ) and f% = 2v\t\ 2 p ail {e d + U). 
Here we assume = and ^ Tf. This reflects 
the fact that the density of states of the left lead is spin- 
degenerate while the right one is spin-split. Assuming a 
constant density of states and a constant tunneling pa- 
rameter t, we have rj? = T^. With this assumption the 
terms with in Eq. ([2]) cancel out,— and one simply 
finds 



n a) -T\ Q<T n a + Tl llT ng\. (7) 



To calculate the current via Eq. ([7]) we need to find 

a fr( 

d 



dt 



where 



r i<r[l - n a - rig. + n n ] - T Wa [n a - n n ] 
+ IW [n a - n-fj.] - T m<7 n u , (8) 
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When = fj, Eq. © becomes 

= r i CT [l - n a - rig] - T WlJ n a + f m<J n s , (13) 

where the terms with n-j-j, cancel out. 

Stationary regime. In this regime (dn a /dt = 0) Eq. 
([T3l) reduces to 



£qict + (r 01cr — r 01o -)ng- 



(14) 



which can be solved for each spin component, thus re- 
sulting in 



roio-rioff + Toi^Toio- 



(15) 



where U a = (T 01cr + T 10lJ )(r 01g . + T ws ) - (T 01cr - 

r i CT )(r ia -Toia)- Using Eq. (I5H into Eq. © we 
obtain 
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(16) 



From Eq. (|16p we can readily evaluate the current po- 
larization £>=(/*— + ij). Next (Sec. Ill) we 
provide some simple analytical results valid when double- 
occupancy is energetically forbidden. Numerical results 
are presented in Sec. IV. 



III. REGIME OF SINGLY OCCUPIED DOT 

As we shall see, the most interesting behavior takes 
place when the channel is completely within the con- 
duction window and e c i + U is far above the Fermi energy 
of the emitter. With this channel configuration we ap- 
proximate f n = 0, and /x, = 1, fa = for eV > and 
Jl = 0, fn = 1 for eV < 0. Using this into the occupa- 
tion and current equations, Eqs. (fT5|) and (fTo]) . we find 
analytical expressions for the first plateau that appears 
in the current and its polarization for both positive and 
negative bias voltage. Equation (| 1 5 j) then becomes 



r»?r T _' 

1 rr 1 /t 



(17) 



where r\ = L, t] = R for eV > and rj = R, r) = L for 
eV < 0. The current of the left lead then becomes 
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(18) 



where ij = R and the + sign corresponds to eV > 0, while 
i] = L and the — sign to eV < 0. The right side current is 
simply given by I R = —1% for a spin-conserving station- 
ary regime. Equation (fTS]) gives the (bias-independent) 
current in the regime addressed here. For the particular 
case of spin- independent tunneling rates, i.e., = T L 
and T R = r fl , we obtain 



I L = If 



2e(r L r fl )/(2r L 
-2e(r L r fl )/(r L 



-v R ). 

2T R ) 



with 



eV > 
eV < 0, 
(19) 

results already known in the 



in accordance 
literature 

Using Eq. (JTHJ) into the definition p = (If — Ij.)/ (If 
L), we obtain the current polarization plateau 
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where r\ = L for eU > and ?y = R for eU < 0. We model 
the tunneling rates by = r[ = T and Ffn\ = r (l± 
p), where p £ [0, 1] is the spin polarization degree of the 



ferromagnetic right lead 2 ^ and To the lead-dot coupling. 
Within this model Eq. (|2"0"|) gives 



P 



0, 

P, 



eV > 
eV < 0. 



(21) 



Thus, when only the level td is within the conduction 
window, the current becomes unpolarized for positive 
bias, while spin-polarized for negative bias. Therefore the 
NM-QD-FM junction functions as a current-polarization 
diode. 



IV. RESULTS 



A. Parameters 



We assume that the dot 
voltage according to = 
counts for asymmetric voltag 
right tunnel barriers i 23 ' 24 e ga t 
voltages. For the numerics 
Ml = 0, [ir = -eV, k B T = 
T = lO/^eV.ii^ In Sees. B, 
metric potential drop across 
In Sec. E we briefly discuss 
x ± 0.5. 



level depends on the bias 
£gate — xeV , where x ac- 
;e drops along the left and 
can be controlled via gate 
we take e ga te = 0.5 meV, 
212 fieV, U = 3 meV and 
C and D we assume a sym- 
the system with x = 0.5. 
the asymmetric case with 



B. Current polarization 

Figure [2] shows the current polarization as a func- 
tion of the external bias eV. We observe that for pos- 
itive bias the current polarization decreases for increas- 
ing bias, reaching zero around eV =4 meV. Conversely, 
for the negative biases we obtain a maximum polariza- 
tion p around eV = —4 meV, confirming the analyt- 
ical result found in Sec. Ill, Eq. ([2"T]) . The voltage 
range for this behavior scales with the parameter U. For 
high enough bias voltages (\eV\ > 7 meV) the polariza- 
tion reaches the same nonzero plateaus for both positive 
and negative voltages. Both the suppression (eV > 0) 
and the enhancement [eV < 0) of the current polar- 
ization are due to the interplay of Coulomb interaction 
and spin accumulation in the quantum dot. Quite in- 
terestingly this interplay affects p differently with the 
bias sign, namely, for direct bias it suppresses p while 
for reverse bias it enhances p.— The suppression of p 
for positive bias results in the zero polarization seen for 
all p values except p = 1. In the half-metallic case 
(p = 1), there is only spin up current flowing in the 
system (/? ^ 0, I? = 0), so the polarization becomes 
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FIG. 2: Current polarization p as a function of the external 
bias. For p values between 0.2 and 0.9, p reaches zero for 
some particular positive bias range, while for the negative 
counterpart it reaches maximum plateaus. 



simply p = (77 - ij 7 )/^ +17) = q/tf = 1. On the 
other hand, for negative bias, the maximum polarization 
plateau changes as p varies. In particular, p attains a 
plateau equal to the polarization degree of the ferromag- 
netic lead, according to Eq. (|2~Tj) . To gain a more detailed 
understanding of the spin-diode effect we investigate next 
the spin accumulation m = rij — and the spin-resolved 
I — V curves as a function of the bias. 



C. Spin accumulation 

Figure [3] shows the spin accumulation to = n-f — nj_ as 
a function of the bias voltage, for distinct polarization 
parameters p. For all the p values considered here we 
note that to < for positive bias and m > for negative 
bias. This spin- imbalance can be understood in terms of 
the tunneling rates between dot and leads. Due to the 
ferromagnetism of the right lead, the rates and be- 
come asymmetric. For example, for p = 0.2 the rates are 
Tf = 12 ^eV, Tf = 8 fieV and Tf = Tf = 10 ^cV. For 

positive bias, becomes the ingoing tunneling rate for 
electrons with spin a and the outgoing tunneling rate. 
Due to the inequality > Tf , the spin up electrons can 
tunnel out the dot faster than they come into it. On the 
other hand, since Tf < Tf, the spin down electrons leave 
the dot slower than they come into it. So on average the 
spin down electrons spend more time in the dot than the 
spin up ones for eV > 0, thus n± > n-\ =>• to < 0. A 
similar reasoning applies to the other p values, except for 
p = for which there is no accumulation. For negative 
bias, and are the outgoing tunneling rates while 
and become the ingoing tunneling rates. As a 
consequence of this interchange, the spin accumulation 



FIG. 3: Spin accumulation m = n-\ — rij, as a function of the 
external bias. For p = (unpolarized lead) there is no spin 
accumulation in the dot. When p increases the spin accumu- 
lation increases as well becoming negative for eV > and 
positive for eV < 0. 

inverts its sign (m > 0). For small p values the spin ac- 
cumulation is essentially an odd function of the bias, Fig. 
® 

When p increases, though, the imbalance becomes 
stronger for positive bias. In particular for p = 1, m 
reaches —1 in the positive bias range corresponding to 
single occupancy (e^ + U > Hl), and a constant plateau 
for all negative bias. This happens because no spin-down 
states are available in the right lead for p = 1, so a spin- 
down electron that enters the dot, coming from the left 
side (eV > 0), cannot leave the dot to the right side. 
Hence a spin-up electron cannot hop into the dot when 
e<j + U > so the accumulation becomes completely 
spin-down polarized for positive bias. For high enough 
bias voltages an additional electron with opposite spin 
can jump into the dot (for both positive and negative 
bias), thus resulting in a suppression (in modulus) of m. 

D. Spin- resolved currents 

In Figure [4] we show the spin resolved currents Jf and 
1 1 as a function of the bias voltage for differing polariza- 
tion parameters p. We observe that for positive bias the 
spin up and spin down currents coincide in the plateaus 
indicated by arrows for any p value. This results in the 
zero current-polarization seen in Fig. [2l In the second 
plateau, though, if attains higher values compared to if, 
which enhances p. The strong suppression of if in the 
first plateau (eV > 0) is attributed to the spin imbalance 
to < observed for the corresponding bias range (see 
Fig. [3]). More specifically, since the dot is predominantly 
spin-down occupied for positive bias, the spin-up elec- 
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FIG. 4: Spin-resolved currents against bias voltage. For the 
plateaus indicated by arrows, Jj lies on top of 7|. This gives 
rise to the p — plateau seen in Fig. [2] For big enough 
polarizations (e.g. p > 0.8) a negative-differential resistance 
range arises in the spin up current. 



o 

"cO 
_N 

o 

CL 

CD 
i_ 
i_ 

O 



c 
o 

"crj 

N 

O 
Q_ 

CD 
i_ 
l_ 

O 



1,0 
0,8 
0,6 
0,4 
0,2 
0,0 



H - 1 .u 

p=u.y /^->. 


x=0.2 


P=0.8 / _ 




p=0.6 AJ 




p=0.4 / ' 




p=0.2 /"> 




p=0 







-20-15-10 -5 5 10 15 20 
bias (meV) 



p=1.0 




-20-15-10 -5 5 10 15 20 
bias (meV) 



trons tend to be more blocked than the spin-down ones, 
thus reducing further if and interestingly locking it on 
top of 1 1 . In contrast, for negative bias we have the popu- 
lation inversion m > 0. This gives a stronger suppression 
of 1 1 as compared to /j, which enhances the difference 
between and I±, and consequently p. When the chan- 
nel td + U reaches resonance (eV < —1 meV) both the 
If and Zj_ plateaus attain values somewhat closer to each 
other, thus reducing the current polarization (see Fig. 

In particular for p = 1 the Zj, is zero for any bias voltage 
since there are no spin-down states available in the right 
lead. The 7| increases slightly (for positive bias) while 
the dot is becoming populated. When the population is 
high enough the Coulomb interaction plays a role and the 
spin up current goes down to zero4^ This gives rise to a 
negative differential conductance at the beginning of the 
first plateau for eV > [see Fig. [4] with p = For 
negative bias (and p = 1) J-j- attains one plateau instead 
of two steps as for the other p values. This is expected 
because the spin-down electrons do not participate in the 
transport in this case, so no Coulomb interaction effect 
arises. Note that for p = 1 the system can operate as a 
mesoscopic current diode^ ' 35 i 36 



E. Effects of the bias-drop asymmetry 

Here we consider the effects of an asymmetric bias 
drop, i.e., x 0.5. As Fig. O shows, the asymmetry 



FIG. 5: Current polarization p as a function of the bias volt- 
age for the asymmetry parameters x = 0.2 (upper panel) and 
x — 0.8 (bottom panel). We observe that the p — plateau 
enlarges for x — 0.2 and shrinks for x — 0.8 when compared 
to the x = 0.5 case. In contrast, the maximum plateau p = p 
reduces for x = 0.2 and enlarges for x = 0.8 when compared 
to the x = 0.5 case. 



in the bias drop gives rise to quantitative, but not qual- 
itative changes. For x = 0.2 the current polarization 
p goes to zero much slower with the bias than it does 
for x = 0.5. This is so because the resonance condition 
6^ < hl (eV > 0), which is necessary to have p = 
[see Sec, piip ]. happens for higher bias when x decreases. 
For negative bias the resonance ed < Hr is reached faster 
(i.e., at lower biases as compared with the x = 0.5 case) 
for decreasing x. This, in turn, translates into a steeper 
enhancement of p which then attains a plateau at p = p 
[see Eq. (|2"Tj) ]. In addition, for x = 0.2 the zero current- 
polarization plateau (eV > 0) enlarges while the maxi- 
mum plateau (eV < 0) shrinks compared to the respec- 
tive x = 0.5 widths. For x = 0.8 the resonance td < Hl 
(eV > 0) takes place faster with the bias when compared 
to the x — 0.2 and x — 0.5 cases. This results in the 
steeper suppression of p and the shrinkage of the zero 
current-polarization bias range. For negative biases the 
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resonance condition e d < [Ar for x = 0.8 is more slowly 
attained with the bias as compared to the x = 0.2 and 
x = 0.5 cases. Consequently, the polarization p reaches 
the plateau at p = p for higher bias voltages (in modu- 
lus). 



where n s is the dot occupation number, with a = —a. 
This occupation can be calculated sclf-consistcntly via 



(A4) 



V. CONCLUSION 

We propose a NM-QD-FM system which operates as 
a diode for the current polarization. More specifically, 
when double-occupancy is forbidden in the system, i.e., 
the channel e d + U is far above the chemical potential 
of the emitter lead, the system carries an unpolarized 
current for positive bias and a spin-polarized current for 
negative bias. This effect is a result of the interplay be- 
tween spin accumulation in the dot and the Coulomb 
interaction. Interestingly, for positive biases the spin- 
resolved currents 7| and ij lock onto the same plateau 
for a particular bias range, thus resulting in p = 0. 
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where the correlation function G^ a (e) is given by the 
Keldysh equation 



G< CT (e)=G^(e)£<G« CT (e). 



(A5) 



The advanced Green function G% a {e) is given by 
G a aa (e) = [G^(e)]*, while £< = i[T^f L + rff R }. In or- 
der to consider the same channel configuration adopted 
in Sec. Ill, we assume a large U (e d + U 3> fi v ) so the 
retarded Green function becomes 



Gw(e) 



(l-n s ) 



and the lesser Green function reads 

r < n= ^{T L JL+T R a fR){l-n s f 



(A6) 



(A7) 



Substituting Eq. (|A7|) into Eq. (|A4|) we have 
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APPENDIX A: NONEQUILIBRIUM GREEN 
FUNCTIONS 

In the main text we have formulated the problem via 
master equation formalism. Here we show that Eqs. (|17p 
and (fT8|) . from which our main result Eq. (j2Tj) directly 
follows, can be obtained via the Keldysh formalism. We 
start with the well known equation for the stationary 

j. 46 47 

curren t °i ' 

V=ie J ^im[<ZAe)-(%Ae)]fv(e)+G< a (e)}}, 

(Al) 

where G r aa , G a aa and G< a are the retarded, advanced 
and lesser Green functions, respectively. To calculate 
these we apply the equation of motion technique and use 
the Hartree-Fock approximation to factorize high-order 
correlation functions in the resulting chain of equations^ 
The retarded Green function becomes 

G r aa {e) - _ lM X yr( y (A2) 
flW(e) - S^(e) 

where S^(e) is the non-interacting tunneling self-energy 
given in the wide band approximation by ££.(e) = 
-ir a /2 = -i(T% + T§)/2, and g aa (e) is the dot Green 
function without coupling to leads, 

- e d - U(l - ns) 



{e-e d )(e-e d -UY 



(A3) 



Now assuming that the dot level e d is completely on res- 
onance within the conduction window between [il and 
for positive or negative bias^ we can integrate Eq. 
in order to obtain 



P7 

n a = =r-(l - rig-), 

J- (T 



(A9) 



where r) = L for eV > or 77 = R for eV < 0. Solving 
Eq. (|A9[) for each spin component we find exactly Eq. 

(ED- 

To obtain Eq. (fl8|) from the Green functions, we sub- 
stitute Eqs. (|A6|) and (|A7[) into the current formula (|A1|) . 
which gives 



de T L „T*{l-n s f{f L - f R ) 
2^ {e-e d y + ^fY{l-n s y 



(A10) 



Solving Eq. (|A10[) with the same assumptions adopted 
previously^ we find 



I L a =± e ^(l-ng), 

J- (7 



(All) 



where + and — signs correspond to eV > and eV < 0, 
respectively. Eq. (|A11[) can also be written as Eq. (TT8"|) . 
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